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We investigate the two-particle interferometry for the particle-emitting sources which undergo 
the first-order phase transition from the quark-gluon plasma with a finite baryon chemical potential 
to hadron resonance gas. The effects of source expansion, lifetime, and particle absorption on the 
transverse interferometry radii -R ou t and -Rsidc are examined. We find that the emission durations of 
the particles become large when the system is initially located at the boundary between the mixed 
phase and the quark-gluon plasma. In this case, the difference between the radii R ou t and 7? s idc 
increases with the transverse momentum of the particle pair significantly. The ratio of \J i?o Ut — -Rsidc 
to the transverse velocity of the pair is an observable for the enhancement of the emission duration. 

PACS numbers: 25.75.-q, 25.75.Gz, 25.75.Nq 



I. INTRODUCTION 

One of the important issues of high energy heavy ion 
collisions is to find and quantify the quantum chromody- 
namics (QCD) phase transition between the quark-gluon 
plasma (QGP) at higher energy density and the hadron 
gas at lower energy density. The initial systems produced 
in the heavy ion collisions at the higher energies of the Su- 
per Proton Synchrotron (SPS) and Relativistic Heavy Ion 
Collider (RHIC) and the energy of the Large Hadron Col- 
lider (LHC) have high temperature and near-zero baryon 
chemical potential. Lattice QCD calculations have shown 
that the transition at the vanishing baryon chemical po- 
tential is a crossover However, it is predicted that 
this crossover will become a first-order phase transition 
at intermediate temperatures and high baryon chemical 
potentials [2-7] . Recently, the search for the evidences of 
the first-order phase transition and location of its critical 
end point have attracted special attention, for instance, 
the RHIC and SPS low energy programs [8l-fl2l] and the 
project of the future Facility for Antiproton and Ion Re- 
search (FAIR) at GSI QMlj. 

Two-particle Hanbury-Brown-Twiss (HBT) interfer- 
ometry is a useful tool for detecting the space-time struc- 
ture of particle-emitting sources in high energy heavy 
ion collisions [l7l - l20j . For the first-order phase transition 
there is a mixed phase of the QGP and hadron gas. In 
the absence of pressure gradient, a slow-burning fireball 
is expected when the initial system is at rest in the mixed 
phase, and this may lead to a considerable time-delay of 
the system evolution [21-26]. It is therefore of interest to 
probe the time-delay for the first-order phase transition 
by HBT interferometry. 



In Ref. [271 ] an HBT analysis technique with quantum 
transport of the interfering pair (QTIP) is developed. It 
takes into account the effects of resonance decay and mul- 
tiple scattering of pions in the sources. In Ref. [28|, this 
HBT technique is used to investigate the source radius 
and lifetime for the spherical systems evolving hydrody- 
namically with the first-order phase transition. In this 
study we use the relativistic hydrodynamics in (2+1) di- 
mension to describe the particle-emitting sources which 
undergo the first-order phase transition. We investigate 
the HBT radii i? ou t, -Rsidc, and -Ri ong [H, H(| for the hy- 
drodynamic sources using the HBT interferometry with 
the QTIP technique. The results indicate that the ra- 
tio of AR 
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of the particle pair v KT is sensitive to the emission du- 
ration of the source. It is large when the system is ini- 
tially located at the boundary of the mixed phase and the 
QGP (soft point). As compared to pion HBT interferom- 
etry kaon HBT interferometry may present more clearly 
the source space-time geometry at the emission, because 
kaons (for instance K + ) can escape easily from the sys- 
tem after their production. By comparing the results of 
the two-pion and two-kaon HBT analyses, we find that 
the particle absorptions and the large expansion veloc- 
ities of the sources after hadronization may change the 
pion HBT radii as functions of the transverse momentum 
of the particle pair. However, the large values of the ratio 
Ai? os /t> KT for the soft point of the first-order phase tran- 
sition can be observed in both of two-pion and two-kaon 
HBT measurements. 

The paper is organized as follows. In section II we 
present briefly the description for the relativistic hydro- 
dynamics in cylindrical coordinate frame. We describe 
the model of the equation of state (EOS) of first-order 
phase transition used in our calculations. The adiabatic 
cooling paths and the space-time evolution of the systems 
are also discussed in this section. In section III we per- 
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form the two-pion and two-kaon HBT analyses, with the 
QTIP technique, for the hydrodynamic particle-emitting 
sources for the initial conditions of the QGP and the soft 
point of the first-order phase transition. The effects of 
source expansion, lifetime, and particle absorptions on 
transverse HBT radii are investigated. On the basis of 
the investigations, we introduce an observable to probe 
the long lifetime of the source for the initial conditions 
of the soft point. Finally, the summary and conclusions 
are presented in section IV. 

II. HYDRODYNAMICAL EVOLUTION WITH 
FIRST-ORDER PHASE TRANSITION 

A. Relativistic hydrodynamic equations in 
cylindrical frame 

The dynamics of ideal fluid in high energy heavy ion 
collisions is defined by the local conservations of energy- 
momentum and net charges (3lL |32[. The continuity 
equations of the conservations of energy-momentum, net 
baryon number, and entropy are 

d^(x)=0, (1) 



8 t N s + d p (N s v p ) + d z (N s v z ) = —N s , (9) 

where E = T 00 , M p = T 0p = T p0 , M z = T 0z = T z0 , 
N b = j% = ?i67, and N s = j° = aj. 

B. Equation of state 

In the equations of motion <[5]) — d9j), there are e, p, v p , 
v z , rib, an d s six unknown functions. In order to obtain 
the solution of the equations of motion, we need an equa- 
tion of state (EOS), p(e, rib, s), which gives a relation for 
p, e, rib, and s. In our model the QGP phase is described 
by a perfect gas of gluons, u, d, s quarks, and antiquarks, 
with the constant vacuum energy B associated with QCD 
confinement [Hj]. The pressure, energy density, and the 
conserved charge density in the QGP phase are given by 

p«=5^ft(T,/i i )-B, (10) 



£«=^ £i (T, Mi )+fl, (11) 



(2) 



Ai m(T,fj,i) . 



(12) 



(3) 



where x is the space-time coordinate of a thermalized 
fluid element in the source center-of-mass frame, T pu (x) 
is the energy momentum tensor of the element, j b (x) = 
n b (x)u p and j p ( x ) = s(x)u M are the four-current-density 
of baryon and entropy (rib and s are the baryon den- 
sity and entropy density), and u p = 7(1, v) is the four- 
velocity of the fluid element. The energy momentum ten- 
sor T^ix) is given by HH, HI] 

T^{x) = [e(x) +p{x)]u> 1 {x)u v (x) - p{x)g» v , (4) 

where p and e are the pressure and energy density of the 
fluid element, and is the metric tensor. 

In the cylindrical coordinate (t, p, <fi, z) frame, g pu = 
diag(l,-l,-p -2 ,-l). The conservation Eqs. CD - © 
can be expressed as 



d t E + d p [(E+p)v p ]+d z [(E+p)v z ]= (E + p), (5) 

P 



d t M» + d p (M p v p + P ) + d z (M p v z ) = M p , (6) 

P 

d t M z +d p {M z v») + d z (M z v z +p) = -—M z , (7) 



d t N b + d p (N b vP) + d z (N b v z ) = N b 



(8) 



where pj(T,//j), e,(T, fn), and rii(T, /x^) are the pressure, 
energy density, and number density of particle species 
i in the perfect gas with temperature T and chemical 
potential {/ii}, Ai is the conserved charge number of the 
particle species i. In our calculations we use the quark 
masses m u = rrid = 5 MeV. m s = 150 MeV and the bag 
constant B = (235 MeV) 4 [33j|. 

For the hadronic phase we adopt the excluded volume 
model [33W35I and consider the particles n, K, N, A, S, 
A, and their antiparticles in the model. The pressure, 
energy density, and the conserved charge density in the 
hadronic phase are given by [33M35I ] 



P 



1 + V J2i n i( T , h) 



h z2iAirii(T,lii 



1 + V J2i n i{T, h) 



where 



[Li = m -V p 



H 



(13) 



(14) 



(15) 



(16) 



Vq = (l/2)(47r/3)(2a) 3 is the excluded volume which is 
assumed to be the same for all hadrons with a = 0.5 fm 
11. 
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For the first-order phase transition, there are Gibbs 
relationships in the mixed phase of the QGP and hadron 
gas. We have T Q = T H , fi N<A = 3p u , ^a jS = 2fx u + fi s , 

/V+,ir Q ,ji- = 0, Hk+,K° = Mu _ Ms; — ) and 

p M = pQ(T; = ^ 5 (17) 



e M = a ^(T, /;„,//,) + (!- a) e H (T,/i„,/i s ), (18) 



n A = anj(T,n u ,n s ) + (1 - a) 7i A (T,fi u ,fj, s ) , (19) 

where /Lt„ and (i s are the chemical potentials of u and s 
quarks, and a = Vq /V is the fraction of the volume oc- 
cupied by the plasma phase. The boundaries of the co- 
existence region are found by putting a = (the hadron 
phase boundary) and a = 1 (the plasma boundary). 

Using the thermodynamical relations of mixed gas one 
can get the entropy densities s and other thermodynam- 
ical quantities, in the QGP, hadronic, and mixed phases 
from Eqs. [flTOJ - (QDL [G3J - CD3], and [(HHISM, 
and get numerically the EOS with the first-order phase 
transition. 



C. Adiabatic paths 

For perfect fluid, the entropy and baryon number of 
the system are conserved during evolution. So the ratio 
of the densities n& and s, nb/s, is a constant. In the 
calculations we take nb/s = 0.06 which corresponds to 
the incident energy about 30 AGeV (36|. The solid lines 
in Fig. [T]show the adiabatic cooling paths for the system 
evolving with the EOS of the first-order phase transition. 
The dotted line is the transition curve between the QGP 
and hadron gas. The mixed phase is on the transition 
curve from the end point of the QGP branch (point 1) 
to the beginning of the hadronic branch (point 2). The 
non-trivial zigzag shape of the trajectory indicates that 
the system has a re-heating in the mixed phase (35L H3] ■ 
The reason is that at a certain point (T, fi) on the phase- 
transition curve, the number of degrees of freedom, and 
hence the specific entropy, is larger in the QGP phase 
than which in the hadronic phase. The temperature must 
increase during hadronization to conserve both the total 
entropy and baryon number simultaneously [37} • 

In Fig. [5] we show the thermodynamical quantity, p/e, 
as a function of e for the system. The ratio p/e reaches 
the minimum at the boundary between the QGP and 
mixed phase, e = e MQ = 1.83 GeV/fm 3 . It is so called 
the soft point of the first-order phase transition. At the 
boundary between the mixed phase and hadronic gas, 
the ratio reaches its maximum. It is named hadroniza- 
tion point. One can see that the ratio retains the values 
smaller than 0.075 in the e regain 0.6 - 2.1 GeV/fm 3 . 
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FIG. 1: (Color online) The adiabatic paths for the system 
evolving with the EOS of the first-order phase transition. 
The dotted line is the transition curve between the QGP and 
hadron gas. The mixed phase is on the transition curve from 
the end point of the QGP branch (point 1) to the beginning 
of the hadronic branch (point 2). 
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FIG. 2: (Color online) The ratio of the system pressure to 
energy density, p/e. It has the minimum at the soft point 1 
and reaches its maximum at the hadronization point 2. The 
ratio retains the values smaller than 0.075 in the e regain 0.6 
-2.1 GeV/fm 3 



D. System evolution 

Using the Sod's operator splitting and RHLLE method 
[3lL l38l - 440j ] . we can obtain the system evolution by solving 
the hydrodynamical equations with the EOS of 

the first-order phase transition. Because the heavy ion 
collisions are full stopped at the energy considered, we 
assume the system is initially at rest within a cylinder 
in the beam direction (z-direction) with the transverse 
and longitudinal radii po and zq. In Fig. [31 we show the 
two-dimension energy density, e{x, z) = J e(x, y, z)dy, for 
the systems at the time t =0, 6, and 12 fm/c. The left 
and right panels are for the systems which are initially 
located in the QGP phase (T QGP = 180 MeV, e^ GP = 
4.12 GeV/fm 3 , /i§ GP = 990 MeV) and at the soft point 
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FIG. 3: The two-dimension energy density, e(x,z), for the 
systems at the time t =0, 6, and 12 fm/c. The left panels are 
for the system initially in the QGP. The right panels are for 
the system initially at the soft point. 
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FIG. 4: (Color online) The average transverse velocity for the 
systems initially in the QGP (a) and at the soft point (b). 




Because there is larger gradient of pressure on the edge 
of the system, the velocity increase more rapidly around 
p ~ pq. At t = 9 fm/c, the decrease of the velocity 
near the center of the system is due to the blast-wave 
expansion which leads to a void in the center region. For 
ICSP, the velocity retains zero in the center region of the 
system even at a larger time because there is not pressure 
gradient in this case. In our calculations, the initial sizes 
for the system with ICQGP are taken to be po = z =4.0 
fin. The initial sizes for the system with ICSP are taken 

~° / £ MQU/3 



to be po = zq 



4.0 x ( £ ?GP , 



=5.2 fin. 



III. HBT INTERFEROMETRY WITH 
QUANTUM TRANSPORT OF THE 
INTERFERING PAIR 

A. Formulas of correlation function 

The two-particle HBT correlation function C(k\, k 2 ) is 
denned as the ratio of the two-particle momentum dis- 
tribution P(ki,k 2 ) to the product of the single-particle 
momentum distribution P(ki)P(k 2 ), 



p(k 1 )p(k 2 y 



(20) 



Using the quantum probability amplitudes in a path- 
integral formalism [23], P(ki) (fcj = [Ej . kA i = 1,2) 
and P(ki, k 2 ) can be expressed as [2^ l4l| - |44j 

P(ki) = f d 4 x p(x)e- 2Xm ^ XK ^\A(xK)\ 2 , (21) 



P(k 1 ,k 2 ) =J d x 1 d 4 x 2 e ""' ^ e 
xp{xi)p(x 2 )\${xix 2 ] kik 2 )\ 2 , 



(22) 



where p{x) is the four-dimension density of the particle- 
emitting source, A(xn) is the amplitude for producing 
a particle at x with momentum k, e ~ 2Xm <f>s{xn^k) ^ g ^. ne 
absorption factor due to the multiple scattering when the 
particle propagating in the source, and $(x\X 2 ; k\k 2 ) is 
the wave function for the two identical bosons, 



(T SP = 142 MeV, el 



SP 



1.83 GeV/fm 3 



SP 



780 MeV). It can be seen that the energy density for 
the system with the initial conditions of the soft point 
(ICSP) decreases more slowly with time than that for the 
system with the initial conditions of the QGP (ICQGP). 
Because there are not the initial velocity and pressure 
gradient in the mixed phase, the expansion of the system 
is slow. Figure 2](a) and (b) show the average transverse 
velocity, (v p (p)) z = J v p (p, z)dz, for the systems with 
ICQGP and ICSP, respectively. For ICQGP, the velocity 
increases rapidly from zero at the beginning (t < 3 fm/c), 
and still increases with time during 3 < / < 6 fm/c. 



$(xix 2 ; kik 2 
1 



{A(x lKl , fc 1 ) J 4(a; 2 «2 ) k 2 )e lk ^ +lk ^ 



V2 

+{A(xik' 2 , k 2 )A{x 2 K t 1 , k%)e 
A(xK,k) = A{xn)e l ^ t{xK ^ k \ 



ik± -xi-\-ik2 -k-2 



h (23) 



(24) 



where S m f(xn — > k) is a phase arising from the source 
collective expansion, which can be described by a long- 
range density-dependent mean- field [13, |44| • 

In our HBT calculations, the identical kaons (for in- 
stance K + ) are assumed to freeze out directly at the 
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hadronization. So, the absorption factor e _2Im ^ 3 is 1 
and <5 m f = 0. The final identical pions (for instance 7r + ) 
include the primary pions emitted at the hadronization 
and the secondary pions from the "excited-state" particle 
decays during the system evolving in hadronic phase un- 
til to the thermal freeze-out. The four-dimension density 
of the pion source can be expressed as [13, El H(| 



p(x) = n 1T (x)d(t — 



^nj(x), (25) 



where rii(x) and r h are the particle number density and 
the hadronization time in local frame, Dj^.^ is the prod- 
uct of the decay rate in time and the fraction of the de- 
cay. For example, Da^tt = Ta x | and D 7r o, II .o_ ! . Tr + 7r - = 
f r n 7r cr(7r 7r — > 7r + 7r~) x 1, where v r is the relative ve- 
locity of the two colliding pions and the cross section 
a(ir°ir — > 7r + 7r _ ) is equal to the absorption cross section 
of 7T+7T- -> 7r°7r° [13. 

When a pion propagating in the source it will subject 
to multiple scattering with the medium particles in the 
source. The absorption factor due to the multiple scat- 
tering in Eqs. (J3TJ and fl22) can be written as (271 l4~0 - |4i | 



e -2Xm<M*) =exp 



o- abs (7ri) rii(x')jd£(x') 



(26) 

where ^ i means the summation for all medium parti- 
cles except for the test pion along the propagating path 
d£(x'), Cabs^*) is the absorption cross section of the pion 
with the particle species i in the medium, and Xf is the 
freeze-out coordinate. In calculations we only consider 
the dominant absorption processes for the identical pi- 
ons, for example the reactions of ir + Tr~ — > tt°tt° and 
tt + N — > A for 7r + , as in Ref. (27|. The pion freeze-out 
temperature is taken to be 110 MeV, which corresponds 
to the energy density et =45 MeV/c [47|. 

In the HBT analysis, we use the Bertsch-Pratt com- 
ponents of the relative momentum q = |ki — ki2 1 of the 
identical particle pair [2§, 30] , q S idc, 9out, and qi on g as vari- 
ables. The correlation function Cxiqside, 9out, Qlong) are 
constructed from P{k\,k2) and P(ki)P{k2) by summing 
over ki and k2 for the (g s ide, 9out, 9iong) bins in a certain 
K T = ||ki - k 2 | T region. The HBT radii R siAc (K T ), 
Routi^r), and i?iong(^T) are obtained by fitting the cor- 
relation functions with the parametrized formula 

C K (q s ide, gout, 3ion S ) = 1 + Xexp[-q 2 ide R 2 ide (K T ) 

-q 2 out R 2 out (K T ) - qf oag Rl ng (K T )} , (27) 

in the longitudinal comoving system (LCMS). Here A is 
called the chaotic parameter. 



B. Results for hydrodynamic sources 

In Fig. [5] we show the two-pion and two-kaon HBT re- 
sults for the hydrodynamic sources for ICQGP and ICSP. 
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FIG. 5: (Color online) The two-pion and two-kaon HBT re- 
sults for the hydrodynamic sources for ICQGP and ICSP. 



It can be seen that there is much difference for the two- 
pion HBT radius i? ut as functions of Kt for the two 
kinds of sources. One decreases with Kt, and another 
almost increase with Kt- When the system is initially 
located at the soft point (ICSP case), the results of i? ou t 
are much larger than those of i? s ido at larger Kt, and the 
ratio i? ut /-Rsidc increases with K t significantly. As com- 
pared to the pion HBT radii the kaon HBT radii exhibit 
more moderate changes with Kt- 

Figure [6] (a) and (b) show the transverse velocities of 
the pion- and kaon-emitting sources as functions of the 
pair transverse momenta Kt- Figure El (c) and (d) show 
the standard deviations of time, a t = y( (t — (t)) 2 ), of 
the particle-emitting sources. One can see that the trans- 
verse velocities of the pion and kaon sources are smaller 
for the system initially at the soft point (ICSP) than 
those for the system initially in the QGP (ICQGP). The 
standard deviations of time enhance very much for the 
sources for ICSP. 

In HBT interferometry, the source HBT radii are re- 
lated to the enhancements of the correlation functions 
at small relative momenta. For an evolving source, 
the source expansion leads to a correlation between the 
particle-emitting coordinate and momentum. It may de- 
crease the transverse emission region for the particle pairs 
with small relative momenta and large Kt- This effect is 
more important in the direction of the transverse momen- 
tum of the pair (out direction), which is boosted by the 
source expansion. Additionally, the source opacity, due 
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FIG. 6: (Color online) Left panels: the transverse velocities of 
the particle-emitting sources as functions of Kt- Right pan- 
els: the standard deviations of time of the particles-emitting 
sources. 



to the absorptions for the particles propagating through 
the center of source (in which the temperatures are higher 
than the hadronization temperature) and by the multiple 
scattering among the particles in the source, may lead to 
a shell emission. This will increase the effect of the de- 
crease of emission region for expanding sources. In Fig 
we show the distributions of the source coordinates pro- 
jected on the transverse out-side plane, for the particles 
with the smaller pair momenta Kt < 300 MeV/c (left 
panels) and the larger pair momentum Kt > 300 MeV/c 
(right panels). The upper four panels are for the system 
initially in the QGP (ICQGP). The lower four panels 
are for the system initially at the soft point (ICSP). For 
Kt > 300 MeV/c, the distributions of the source coor- 
dinates are more concentrated in r Q > regions. For 
Kt < 300 MeV/c, the annular distributions for kaon in- 
dicate that the sources are almost transparent for the 
kaons emitted later. We will see it is that the source ex- 
pansion, lifetime ( ~ at), and particle absorptions lead to 
the differences of the transverse HBT radii i? ut for the 
two kinds of sources for ICQGP and ICSP. 



C. The effects of source expansion and lifetime on 
transverse HBT radii 

In HBT interferometry, the difference of the transverse 
HBT radii in out and side directions includes the impor- 
tant information on the source expansion and lifetime. 



m de is given by [1J, |48|, |49j 



R„ 



R. 



side 



= ((r 



If) 



>1) 



(28) 



where ( • • • ) denotes the average for the space-time co- 
ordinates of the source, f a and f s are the biases of the 
source spatial coordinates related to their average values 
in the out and side directions, t is the bias of the source 
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FIG. 7: The distributions of the source coordinates projected 
on the transverse out-side plane. The left and right panels are 
for the smaller and larger particle pair momenta Kt < 300 
MeV/c and Kt > 300 MeV/c. The upper four panels are for 
ICQGP. The lower four panels are for ICSP. 



and 



is the 



time coordinate related to its average, 
transverse velocity of the particle pair. 

In order to examine the effects of source expansion 
and lifetime on the transverse HBT radii i? out and i? s id , 
we investigate next the two-pion interferometry for the 
simple sources with a constant temperature 100 MeV and 
the Gaussian space-time distributions as 



dN 
d 3 rdt 



cx exp 



y 



2Rl 



2r 2 



Ri < Vx* 



y 



z 2 < R 2 



(29) 



We take Rt = Rl = 5 fm, R% = 10 fm, and assume that 
the sources have the radial velocity 



V0 R~ 2 - 



(30) 



Here r, R±, and vq are three free parameters. We taken 
r = 6 and 12 fm/c for the sources with shorter and longer 
lifetimes. For a shell source R\ is taken to be 5 fm. For 
static and expanding sources, vq is taken to be and 0.8 
respectively. Because there are not correlations between 
spatial coordinates and time for these sources, Eq. 
reduces to 



p2 r>2 

^out "* 



side 



(0-<r-)+cr t X T , 



(31) 



7 



where of ={(t- (t)) 2 ) = (tt - 2)t 2 /tt a 0.363r 2 . 

Because of source expansion and opacity the difference 
of the variances in out and side directions, (ff) — (f 2 ), 
is not zero even for the source with transverse symme- 
try. It is negative and decrease with Kt- On the other 
hand, the right third term in Eq. (f3"Tj). erf V ^ T , is positive. 
It increases with Kt and becomes important when the 
source lifetime r increases. 

In Fig. [51 we show the transvers e HBT radii i? ou t and 
i? s ido and Ai? os = yj R 2 ut — R 2 ide for the sources with 
r = 6 and 12 fm. The symbols o, v, and A are for 
the static Gaussian source (v = 0, R± = fm), ex- 
panding Gaussian source (vq = 0.8, Ri = fm), and 
expanding shell source (vq = 0.8, R\ = 5 fm). The 
dashed lines in the bottom panels are the results of a t v KT 
{v KT = K T /E K , E K = {E 1 +E 2 )/2.) In Fig. M we show 
the distributions of the source coordinates projected on 
t — r s plane. The panels (a), (b), and (c) are for the 
static Gaussian source, expanding Gaussian source, and 
expanding shell source for the smaller pion pair momen- 
tum K T < 300 MeV/c. The panels (a'), (b'), and (c') 
are for the static Gaussian source, expanding Gaussian 
source, and expanding shell source for Kt > 300 MeV/c. 
For the static sources, the results of R s ide are almost a 
constant and R ou t increases with Kt- Because there is 
not the effect of source expansion, (rf) — (r 2 ), and the 
results of Ai? os are consistent with those of atV KT . For 
the expanding and shell-emitting sources, the source ex- 
pansion and shell emission change the distributions of 
the source coordinates. It leads to the decreases of i? s idc 
with Kt- Although i? ut increases with Kt in the small 
Kt region for the sources with larger lifetime r, this in- 
crease will be counteracted at large Kt by the effects 
of the source expansion and shell emission. In these 
cases, the results of Ai? os are smaller than the values 
of a t v KT at larger Kt- From Fig. [§] one can see directly 
that the coordinate distributions of the static sources for 
the smaller and larger pion pair momenta are almost the 
same. However, the coordinate distributions of the ex- 
panding sources for Kt > 300 MeV/c are more concen- 
trated in r Q > regions as compared to the corresponding 
distributions for K T < 300 MeV/c. 

For hydrodynamic sources, there are also correlations 
between source spatial coordinates and time. We will 
discuss the effect of the correlation between r a and t on 
Ai? os in next subsection. 



D. Characteristic quantity for long source lifetime 
for ICSP 



Because of the correlation between source spatial co- 
ordinate r a and time for hydrodynamic sources, Eq. 
becomes 



R n 



R%de = ( rl) - ( f 2 s ) + afv 2 KT -2(rJ)v KT . (32) 



For positive or negative (f t), the right last term in Eq. 
(|3"2"]l will decrease or increase Ai? os with Kt increase. 
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FIG. 10: The distributions of the space-time coordinates of 
source points projected on r — t plane for ICQGP ((a) and 
(c)) and ICSP ((b) and (d)). The dashed lines are for the 
average values of r and t. 



In Fig. [TU]we show the distributions of the space-time 
coordinates of source points projected on r Q — t plane for 
ICQGP (panels (a) and (c)) and ICSP (panels (b) and 
(d)). The dashed lines are for the average values of r Q 
and t, which divide the plane into four regions. In regions 
I and IV, (f t) < 0. In regions II and III, (f a t) > 0. 
For pion, because of particle absorption the distributions 
in region III are less than those in region I. So the values 
of r t averaging over all r D < (r a ) regions (I and III) are 
negative. In r a > ( r a ) regions, the distribution for the 
pion source for ICQGP (panel (a)) is much different from 
that for ICSP (panel (b)) because of the larger transverse 
expansion of source for ICQGP. In region II of the panel 
(a), the wider distribution for the pion source for IC- 
QGP leads to a greater contribution to ( f 1 ) . So the 
value of f t averaging over all r > (r a ) regions (II and 
IV) is positive for ICQGP. For ICSP, the distribution 
for pion source in region II is less than that in region 
IV. The value of f Q t averaging over all r Q > (r a ) re- 
gions is also negative for ICSP. The calculations indicate 
that for pion (f t) =0.16 and -8.30 fm 2 /c for ICQGP 
and ICSP, which are consistent with the above discus- 
sions. For kaon, because of the high transparency and 
low transverse expansion of sources, the distributions are 
approximately symmetric about ( t ) . The values of ( f Q t ) 
are small for the kaon sources. 

In Fig. [TTJ we show the results of Ai? os (symbols o) 
and the products atv KT (symbols *) of the average emis- 
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FIG. 11: (Color online) The results of AR os (o) and <7tD K 
(*) for the hydrodynamic sources for ICQGP and ICSP. 



sion durations at of the particles and the pair transverse 
velocities u KT for the hydrodynamic sources for ICQGP 
and ICSP. It can be seen that except for the pion results 
in panel (a), the results of the Ai? os and a t v KT are almost 
consistent. 

Inspired by the consistences of the results of Ai? os and 
<J t v KT for ICSP, we introduce the quantity 



Ai?o 



out 



7->2 

^sidc 



K T /E. 



K 



(33) 



to describe the character of the long lifetime of the 
sources for ICSP. It is an experimental observable. 

In Fig. [T2] we show the results of a for pion and kaon 
for the hydrodynamic sources for ICQGP and ICSP. The 
larger values of a for the soft point of the first-order phase 
transition are observed in both of the pion and kaon in- 
terferometry measurements. The a values for the pion 
source for ICQGP are much smaller than the average 
value at =8.4 fm/c at large Kt, because of the large 
transverse velocities of the source and the positive values 
of ( f a t ) . In this case a cannot reflect the real lifetime 
of the source. At small Kt-, the larger values of a for 
the pion source for ICSP are due to the large negative 
values of (f t) as well as the small transverse velocities 
of the source in this case. The errors of b exhibited in 
Fig. [12] are only from the statistic errors of i? ut and 
-Rsidc related to the HBT parametrized fits. In high en- 
ergy heavy ion collisions, there are other effects which 
may bring uncertainty to the observable, for example the 
non-equilibrium dynamics during the decay of resonances 
after the hadronization. Further investigations on these 
effects will be of interest. 
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undergo the first-order phase transition from the quark- 
gluon plasma with finite baryon chemical potentials to 
hadron resonance gas. The effects of source expansion, 
lifetime, and particle absorption on the HBT radii are 
examined. For pion, the large transverse expansion of 
the source for ICQGP decreases the HBT radii i? ut and 
R s ide at large transverse momentum of particle pair Kt- 
The source has a long lifetime and small expansion when 
the system is initially located at the boundary between 
the mixed phase and the QGP (soft point). In this case, 
the difference between the transverse HBT radii i? ou t 
and i? s idc increases with significantly. The ratio of 
\J i?Q Ut — Rg idc to the transverse velocity of the particle 
pair u KT , tr, is an observable for probing the long lifetime 
of the source for the soft point of the first-order phase 
transition. As compared to pion HBT interferometry 
kaon HBT interferometry may present more clearly 
the source space-time geometry at the emission. The 
larger values of a for the soft point of the first-order 
phase transition can be observed both by two-pion 
and two-kaon HBT measurements. Further investiga- 
tions on other effects on the observable will be of interest. 
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